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Local field potentials (LFPs), measured using extracellu- 
lar electrodes, provide information about local neuronal 
population activity. The development of multi-electrode 
arrays (MEAs) that measure LFPs simultaneously at dif- 
ferent network locations has led to a renewed interest in 
the LFP as a measure of neuronal network activity [1]. 
While there is a well-understood relationship between 
neuronal current-sources and LFPs, the complexities of 
current generation in spatially organised neuronal net- 
works mean that relating LFP measurements to the 
underlying network dynamics is non-trivial [1]. Forward 
modelling studies have shown that the use of multicom- 
partment models incorporating dendritic structure is 
crucial for simulating spatially realistic LFPs [2]. How- 
ever, most forward models of neural network dynamics 
use simplified single-compartment neurons, because 
they are easier to parameterise and much less computa- 
tionally expensive to simulate. 

Here, we describe two contributions to the effort to 
understand LFP generation in neuronal networks. First, 
we investigated reduced compartmental models to assess 
their suitability for producing reahstic spatial LFP charac- 
teristics. We show that the chosen reduction method [3] 
creates compartmental models that preserve essential spa- 
tial LFP features, while being computationally cheap 
enough for large network simulations. Second, we devel- 
oped the Virtual Electrode Recording Tool for Extracellu- 
lar potentials (VERTEX): a tool for simulating LFPs in 
spatially organised, spiking neural networks. We imple- 
mented VERTEX in Matlab (Mathworks Inc., Natick, MA, 
USA) as it is widely used throughout the neuroscience 
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community, and we provide a simple interface based on 
parameter specification so that VERTEX can be used even 
by non-programmers. Code vectorisation and parallelisa- 
tion methods allow VERTEX to simulate networks con- 
taining hundreds of thousands of reduced compartmental 
neurons in reasonable time. We demonstrate how VER- 
TEX can be used to simulate LFPs at hundreds of loca- 
tions in a large-scale, spatially organised cortical network 
model, and how these results can be used in conjunction 
with experimental data, by comparing simulation results 
with MEA recordings of persistent gamma oscillations 
made in macaque neocortical tissue in vitro. 
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